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Abstract 

A new test apparatus is reconstructed and is applied to investigate static and 

UrlZmtrT?eT?f° f lnnUllr 8C1 ' S ' elked by ' W ° Ph “ e "° W <*" *" d “<■"«) for 

urbopumps. The fluid forces acting on the seals are measured for various parameters 

“ V “ d rlt, °' ' he P' eswil > velocity, the pressure difference between the inlet and 
on.Ie. of the seal, the whirling amplitude, and the ratio of whirling speed to spLtag 

Metis iniLTgaTed ft " e r °' PMlme ‘ eIS 0 " ‘ he SllUc “ d <*■“* character 

Phte t j as th d eXpCnm ™ W IC8 " hs - As 1 ""l*. with regard to the two 
phase flow, as the void ratio increases, the flow induced force decreases. Another dv 

nanttc characler.st.es of two phase flow is as almost similer as that of the mono phase 


1. Introduction 

sure ,-d r b T aChineS SUCh ^ a PUmP tends t0 be ° perated at hi 8 h s P eed ^d high pres- 
sircnmstances theflTd T* ° f performance requirements. In the special 

seal so I" CaUSeS “ inStabmty ° f thC r ° tating machine beCause * 
oncontacting seal. So, it is necessary to make the dynamic behavior of two phase flow 

a^d still m° SUPP y aCCUIate data in order to P redict a » d prevent the unstable vibration 
and still more necessary to design a stable rotating machine 

Many theoretical and experimental analyses on the characteristics of noncontacting 

sea s Childs e^al (L T ^ A * “ expeiimei *tal research on annular 

eals, Childs t al. (Reis. 2,3) mvest.gated fluid forces acting on the seals and their char- 
acteristics with a test apparatus which has an eccentric rotor. Nordmann (Ref 4) used 

“ ,mp : lse t0 '“ *PP«*‘» •« carry out an experiment on .he c aracte isto 

o annular seals, Kaneko e. ah (Ref. 7), Kanki e. al. (Ref. and HoH e. al M 

) also reported Iheir test results on the seals. All these results are important for the 
research on monophase flow annular seals. important lor the 

The authors also investigated for the many types of seals, that is annular seal 
parallel groove seal, spiral seal and tri angle hale seal and also long annular seal for inves’ 
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Of the parameters such as whirling amplitude, pressure difference between the in et and 
outlet of the seal, preswirl velocity, rotating and whirling speeds, and their directions is 
investigated. Then the effect of two phase flow on the dynamic characteristics of annular 

parallel seal is investigated. 


2. Nomenclature 


F X , Fy 
F r ,F t 

R XX > Fyy, K x y, liyx 
ClI) Cyy,C X y,Cy X 

Af.rH Afyy, Af T V ■ Afyi 

r,t 


p{4>,z ) 

P(4>) 


Pes 

A P 


u 

n 


<t> 

R,D 

L 

C 

Vt 

V ti 

V a 


Q g , Q i 

P ) Pi atm 


: Fluid forces in x and y directions 
: Fluid forces in r and t directions 
: Stiffness coefficients 
: Damping coefficients 
: Inertia coefficient 

: Fixed coordinates (z: axial direction) 

; Radial and tangential coordinates 
: Pressure distribution 
: Average pressure in axial direction 
: Inlet pressure to the seal 
: Outlet pressure from the seal 

: Pressure difference between inlet and outlet of the seal 
: Spinning angular velocity of the rotoT 
: Whirling angular velocity of the rotor 
: Whirling eccentricity of the rotor 

: Phase difference between principle force and displacement 
: Rotor radius and diameter, respectively 
: Seal length 
: Seal clearance 
: Preswirl velocity 

: Preswirl velocity without a rotating motion 
: Fluid average velocity in axial direction 
: mean void ratio 

: flow rate of gas and liquid, respectively 
: average pressure in the seal and outlet pressure 


3. Experimental apparatus 


3.1 Test seal apparatus 

Figs 1 and 2 show assembly of the test apparatus and the layout of the test facility, 
respectively In Fig.l, a working fluid, that is, water and gas, is injected through three 
pXTswirl passages to accomplish the different inlet swirl velocities shown in cross 
section R. The water passes through the clearances of the seal and flows to ou e s in 
both sides of the housing. Cool water and gas is continuously supplied to the tanks in 
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order to maintain constant temperature. The inlet part consists of four tubes connected 
will, swirl passages for injecting the water and gas, and swirl speed is adjusted by a 
valve for air and two valves for water in order to obtain the arbitrary swirl velocity in 
e range of 0 - 15.0m/s. These valves are used in order to change the void ratio. 

he seal assembly consists of a seal stator and a sea] rotor. The seal rotor of 
w ic lameter is 70mm, is connected to the motor by a flexible coupling. The motor 
is controlled in the speed range 0 - 3500r.p.m. by an electric inverter, which also selects 
he rotational direction. The seal stator has holes to measure the dynamic pressure in 
the seal, as shown in cross section C. Both long and short seals have three holes in axial 
irec ion an our holes in x and y directions. The dynamic pressure is measured by 
the strain gauge type pressure gauge shown in cross section C. The fluid dynamic force 
acting on the stator is also directly measured by the load cells shown in cross section C 
for comparison with the data of the pressure transducers. 

The bearing assembly has two ball bearings to make spinning and whirling motions, 
o make a whirling motion, an inside sleeve and an outside sleeve, which have a 0.05mm 
eccentricity to each other, are attached between the two bearings as shown in cross 
section A of Fig.l. The two eccentric sleeves can be rotated, relatively. So an arbitrary 
eccentricity can be adjusted in the range of 0 - 0.1mm. The sleeves of both sides are 
driven by a motor through the timing belts. The motor can also be controlled by an 

e ectnc inverter in the rotating speed range of 0 - 3500r.p.m., and rotational direction 
can be selected by an inverter. 

The two phase flow is made as shown in Fig.3, where gas is mixed by higher 
pressure a liquid. The inlet pressure of two phase flow is measured just before the seal 

The void ratio in two phase flow can be changed the pressures of water and gas by 
of valves and regulator. 6 J 


3.2 Measuring instruments and methods 

The measuring procedure illustrated in Fig.4 consists of four kinds of physical vari- 
a les: that is, rotating and whirling speed, dynamic pressure in the seal, seal forces and 
displacement of the rotor. Signals from measuring instruments are recorded by a data 
recorder and analyzed by a computer. 

The rotating and whirling speeds are measured by eddycurrent type pulse sensors 
and digital counters. 

Dynamic fluid force can be directly measured by a load cell as shown in Fig.5. The force 

measure y the load cell is calibrated to revise the influence of the O - rings and the 
inertia of the seal. 6 

An eddycurrent type displacement sensor is used to measure the displacement and 
this displacement is used as a reference signal to obtain the phase difference between the 
displacement and the flow induced force. 

A in , TI ' e T < T“ rc Iecorded to » <!*»» recorder and seal to the computer through au 
A/D converter. In the computer the pressure values are integrated in the circumferential 
direction to obtain the fluid force, then the characteristic coefficients are calculated 

shown in Re « Ti° ‘ c "* i? ‘° mei “ iUIe lhe preswirl velocil J in lhe ^ i-Iet, as 

instead of the usual' U * he 

Ariel flow velocity in seal is determined from the outlet flow velocity directly mea- 
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suied by a pitot tube. 

The flux of water and gas is obtaind by measuring the leakage of unit period by 
cylinder. 

The void ratio in the two phase flow is decided as average value between inlet and 
outlet of seal. Here, it is supposed that air is ideal gas and slip ratio between water and 
gas is 1 because the seal is set horizontally . At this time, we obtain flowing equation 
from the flux of water and gas and the pressure which averages one of inlet and outlet. 


P itrnQg 

r (q, + zyrQg) 


0) 


Seal dimensions that is used for the test and experimental conditions are shown in 
the table 1,2. 


3.3 Calibration and measuring error 

As a preliminary test, the static and the dynamic characteristics of measuring in- 
struments, i.e. the pressure transducer, the load cell and the pitot tube, are calibrated. 

Dynamic calibration of the load cell was done in various working conditions. A 
periodic force was excited on the seal by a shaker, and the force was simultaneously 
measured by strain gauges attached to the shaking rod and by the load cell fixed on the 
other side. These data were used to calibrate the measured data in the real test. 

A pitot tube set was used to measure the inlet preswirl velocity. Its pitot tube 
coefficient for calibration was determined by means of a standard pitot tube. 


4. Calculation of fluid force and characteristic coefficients 

It is assumed that the motion equation of the rotor system is represented as follows 


[m] {x} + [C] {x} + [K] {x} = {/(*)} - {F} 

where [M], [C], [K] are the mass, the damping, and the stiffness matrix and { F(t) 
} is the force vecter. { F } = { F x F y }* is the fluid reaction force acting on a rotor 
and is represented as a linear function of rotor displacement, velocity and acceleration, 
as follows. 



My X 


Ml » l / f \ + f CXI 

Myy \ V ) . Cyx 




( 2 ) 


For a small whirling motion about a center position, the relation of the coefficients 
may be expressed by 
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( 3 ) 


M xx - M yy , M yx = M xy = 0 , C xx = Cyy, 

C V X = ~C X y, K XX = Kyy, J(y X = - K X y 


where ‘he cross-coupled inertia coefficient is neglected because it is negligible small 

force easily TteTlat mi Tr em T"” 8 Wi “ ' he rotor ^“P"* 1 “> >"%« the fluid 
force easily. The relation between a fixed coordinate system and a rotating coordinate 

system ,s illustrated Fig.7, where F, and F, are represented by the radid force F 

and the tangential force F,. Details to obtain the fluid force is shown in (Re 11) so 
these are not written here. vivei.uj, so 


5, Experimental results and discussions 

5,1 Visualization 


We visualize a flow patteern of two phase flow in the part of seal made of acrvle 
As a result, we can observe that bubble of air in two phase flow is fine and two phase 
flow ,s homogeneous bubbly flow because of high rotating speed of the rotor. 


5.2 Relative Uncertainly 

Since the discrepancy of the fluid forces measured by the pressure transducer and 

^nr^ltnln 7d I 0 ' " 16 V ° id ' a,i0 S=0 “ Sh ° W ” in Fie ' 8 ' h is d « ided ‘hat the 
load II ‘ ' es a a is very good. So in this test the fluid force measured by the 

load cell „ used as the experimental results in this discussion 

ratio J-0 9 a 5 „d O 0 «' p"" “n' 0 ' 5 '/ " Uid me “" red '»*-* cell fo, mean void 

°' 45 ' ft0m ,hWe daU “ “ k "°”” ‘ h *‘ P"-e flow has inherent 

the nh "TT* P ' 0CeSS “ lhe VOid rllio w “ increased, random vibration due to 
p use ow ecame very large. So the variance of the measured vibration data 
become large, but these are not the measurement error. 


5.3 Effect of inlet pressure 

F ;V°, Sh T the 6ffeCt ° f inlet pressure on the radial and tangential forces From 

~ ,niei ~ 
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5.4 Effect of spinning speed 

Figs. 11(a) to (d) show the effect of spinning speed on F r /e and F/e for different 
void ratio. Fig.(a) shows for void ratio zero in order to compare the result with another 
void ratio. Figs.(b) to (d) show F r /e and F/e for the void ratio <>=0.25, 0.45 and 0.7 in 
the case of various spinning speeds. From these results it is known that the fluid force 
decreases with the increase of void ratio. The variation of radial force increases with 
increase of void ratio. These phenomena are seemed to be occurred by the compressibil- 
ity of the gas. For the tendency of the tangential and radial forces for rotating speed is 
same, and as rotating speed is increased, the fluid force increases. 

Fig. 12 shows the equivalent spring coefficients and damping coefficients and added 
mass. From these figures it is known that as the void ratio is increased, the spring co- 
efficients and damping coefficients decrease as discussed m the result of F r /e and Ft/e. 
Fig. 13 shows the whirl frequency ratio, 


/ = K iy /c xx n 

and fluid force F/e, which is important to evaluate the instability force ; 


F/r = \\JW+Tt 

From this figure it is known that the whirl frequency ratio decreases at tf=0.25 and 
then increases for increase the void ratio. The whirl frequency ratio has large difference 
at a=0.7, because the variation of force due to two phase flow becomes large and the 
error of amplitude and phase may be large. In this void ratio region the fluid force F/e 
becomes small, then it seemed to be small effect on the rotor system. Conventionaly, the 
pump rotor is supported by bearings and seals in the operational condition, and the seal 
force act as a stabilizing force in the Totor system. In order to demonstrate the effect of 
two phase flow, we assume that a rotor system is stable and in this system the bearings 
act to unstable and the seals act to stable. In this case, if the seal force is decreased by 
the two phase flow as shown in Fig.13, this stabilizing seal force decreases and the rotor 

system may become unstable. 


5.5 Effect of preswirl velocity 

Figs. 14(a) to (c) shows the effect of the whirl ratio v/Sl on F-/e and F/e for various 
preswirl velocities. In this case preswirl velocity ratio is 1.6 for w=2500rpm. The results 
illustrate that the effect of preswirl on the fluid force is same tendency as the case of 
void ratio o=0 but it decreases as the void ratio increases. 

Fig. 15 shows the whirl frequency Tatio versus preswirl velocity. It is known from 
this figure that the effect of preswirl velocity decreases as the void ratio increases. There- 
fore two phase flow of high void ratio does not affect on instability of rotor system. 
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6. Conclusions 


The two phase flow induced force through annular seal are experimentally inves- 
tigated for setting the parameter void ratio, rotating speed, whirl/spin ratio, preswirl 
velocity, inlet pressure. Then the following results are obtained, 

(1) The void ratio is increased, random vibration due to the two phase flow becomes 

very large. 

(2) The fluid force increases with increasing the inlet pressure and its tendency is not 

change with increasing the void ratio. 

(3) The fluid force decreases with the increase of void ratio. These phenomena is 

occurred by the compressibility of the gas. 

(4) Whirl frequency ratio which is proposed by D. Childs decreases a little at a=0.25 

but in practical, it is not so much changed and as the fluid force decreases in high 
void ratio, instability force decreases. 

(5) Theeffect of preswirl velocity decreases as the void ratio increases, so this effect 

dose not act as instability. 

(6) Since the seal force is decreased by the two phase How, this stabilising seal force 

decrease, and rotor system may become unstable. 
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Fig.l Test apparatus assembly 
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Table 1 Seal dimension 



Length ( mm j 


L/D 

Clearance (mm) 

Plain Seal 

70 

70 

1 

a t — 

0.5 


Table 2 Experimental condition 


Rotating speed 

w(rpm) 

500 ~ 3500 

Whirling speed 

H(rpm) 

±600 ~ ±2400 

Whirling amplituide 

e(/j, m) 

50 

Preswirl velocity without a rotating motion 

Vts 

-15 - 15 

Temperature of water 

T(°C) 

19 
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Liquid Gas 



Fig.3 Detail of gas-liquid mixing device 



Fig.5 Detail of force measurement 


Displacement pickup Load cells 



Fig.4 Analysis procedure 



Fig.6 Detail of swirl velocity measurement 
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Fig. 9 Time 










(a) P in = 120kP a 



VhJrJ/SpJn 


(b) P tn = 470kP„ 


Fig. 10 Radial and tangential forces for the constant void ratio 
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Fr/e (kN/mm) b t ^/ e > Fr/e (kN/mm) 
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* 
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* 
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0 
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(c) a = 0.45 


1 a = 0.7 

€ s 0.05 mm 

Vts = 0 m/s 

2 1 u -* 1 1 i 

-2 -J.S -I -.5 0 .5 1 J.s 
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(d) a=0.7 


Fig. 11 Radial and tangential fluid forces for void ratios 
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Fig. 12 Fluid force cofRcients for void ratios 



0 .2 .4 .6 ,0 


void ratio 

Fig. 13 Whril frequency ratio and fluid force versus void ratio 
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ntial forces for preswirl velocity 
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